Magnetic properties of graphene quantum dots 
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Using the tight-binding approximation we calculated the magnetic susceptibility of graphene 
quantum dots (GQDs) of different geometrical shapes and characteristic sizes of 2-10 nm, when 
the magnetic properties are governed by the electron edge states. Two types of edge states can 
be discerned: the zero-energy states (ZES) located exactly at the zero-energy Dirac point, and the 
dispersed edge states (DES) with the energy close, but not exactly equal to zero. DES are respon- 
sible for the temperature independent diamagnetic response, while ZES provide the temperature 
dependent spin Curie paramagnetism. The hexagonal, circular and randomly shaped GQD contain 
mainly DES and, as a result, they are diamagnetic. The edge states of the triangular GQD are 
of ZES type. These dots reveal the crossover between spin paramagnetism, dominating for small 
dots and at low temperatures, and orbital diamagnetism, dominating for large dots and at high 
temperatures. 
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I. INTRODUCTION 

In past years the special attention has been paid to fab- 
rication of graphene quantum dots (GQDs)ii^ susceptible 
to be used for magnetic field-controlled spin-electronic 
logic gates^. However the origin of magnetism in such 
structures still remains unclear in spite of number of re- 
cent studies of orbital^^i^ and spir*i2rJ^ properties. This 
concerns, in particulary, the dependence of magnetic re- 
sponse on the size and on the shape of GQD. 

Landau diamagnetism in perfect infinite graphene 
sheet was first studied by McClurei^"— and more recently 
in Refs. [r5420| . were the singular behavior of suscepti- 
bility was found when the Fermi energy approaches the 
Dirac point at zero temperature. This peculiar behavior 
around zero energy also takes place in the cases when 
disorder-provided band is present for infinite graphene 
and ribbonsi^"— On the other hand, the presence of the 
edge states with energy around zero is a signature of the 
graphene nanoflakes with various terminations and most 
notably with the zig-zag edgesi^"— 

The number and the properties of edge states are sen- 
sitive to the geometry of the GQDi^"— Since diamag- 
netism of graphene occurs due to the electronic states 
with the energy near the Dirac point, it is natural to 
assume that the edge states should make a dominant 
contribution to magnetism of graphene nanoflakes and 
the geometry of GQD will play an important role in the 
diamagnetic response of the nanostructure. 

In this paper we study the hexagonal, circular, triangu- 
lar and random GQDs and identify two types of the edge 
states. First, there are the dispersed edge states (DES) 
whose energies are distributed in the range of 2A around 
the Dirac point, with the value of A being inversely pro- 
portional to the size of the GQD. Secondly, there could be 
highly degenerate exactly-zero-energy states (ZES). The 
DES are appropriate to the hexagonal, circular and ran- 
dom GQDs. Their energies are sensitive to the applied 
field that induces the edge currents. These states pro- 
vides the orbital diamagnetic response of the nanoflakes. 



The number of ZES, that are mostly present in triangu- 
lar GQDs, can be found exactly from the graph theory^. 
Their origin is purely geometric and their location does 
not change as function of the applied magnetic field. 
Therefore ZES do not contribute to the diamagnetism 
of GQDs, but they can be occupied by the electrons 
with unpaired spins and provide the paramagnetism of 
the systen*ii. Studying the edge-state-provided orbital- 
diamagnetic and spin-paramagnetic response of GQD we 
predict the possibility of crossover between paramagnetic 
and diamagnetic response of GQD as a function of their 
shape, size and temperature. 

After this work was completed, we became aware of 
the preprint^ where the similar research was done. The 
main difference between our results is that, we con- 
sider the dots of smaller size and at low temperatures 
ksT <C A, where the edge-states-provided diamagnetic 
peak is broaden b y t he size effects and is temperature 
independent. Ref. [33| mainly addresses the temperature 
effects relevant for the GQDs with bigger sizes and small 
values of A when diamagnetism is mainly of the bulk 
origin. 



II. THE MODEL OF GRAPHENE QUANTUM 
DOTS 



We use the simplest nearest-neighbor tight-binding 
approximation where the properties of conducting tt- 
electrons of graphene are described by the Hamiltonian 



(ij) 



(1) 



where c|, Ci are the creation and annihilation electron 
operators and Si is the on-site energy. In what follows 
we do not consider any on-site disorder and set Si = 
0. The hopping matrix elements between nearest- 
neighbor carbon atoms account for the magnetic field via 
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the Peierls substitution, 

Uj = 70 exp 1^ ^ A • dlj , (2) 

where A = (0, Bx, 0) is the vector potential of magnetic 
field and the zero-field hopping was taken as 70 = 3.0 eV. 

The graphene flakes were selected of hexagonal, circu- 
lar, triangular and random shapes with mostly zig-zag 
edges. The contour of the random shape nanostructures 
has been defined in polar coordinates by: 

r{e) = ro + £ [Ak cos{ke) + Bk sm{ke)] (3) 

k=l 

where rg is the constant average radius that defines the 
typical size of the fiakes, and Bk arc the random 
numbers with amplitude not exceeding ro/3. In order 
to have realistic variation of the flake edge on the scale 
of the lattice constant, the maximum number of harmon- 
ics /cniax has been chosen to be of the order of ro/a, where 
a ~ 2.461 A is the lattice constant. Typically, this num- 
ber was about 25. 

Direct numerical diagonalization of Hamiltonian ([1]) 
gives the fleld-dependent energy levels En{B) and cor- 
responding on-site amplitudes (pn.i of the wave function. 
The orbital energy of the 7r-clcctrons at zero temperature 
as a function of the chemical potential fi and magnetic 
held B is given by 

U{B,fi) = 2 En{B), (4) 

n 

where the factor of 2 is the spin-degeneracy of the levels. 
The low-temperature diamagnetic susceptibility per unit 
area has been calculated as 

1 \ d^UiB,e) ] 

XW = — jr^ — , 5 

^ L dB^ J 

where a = \/3a'^N/4 is the area of graphene flake con- 
taining N carbon atoms. 

III. BULK AND EDGE STATES 

In what follows it will be convenient to distinguish the 
bulk and the edge electronic states using the following 
geometrical criterium. For a given state with the energy 

En we ascribe the intensity In'' ~ J2r <b, IV'n.iP of the 
electronic states located within the circle of radius R to 
the bulk part of the total wave function intensity, whereas 

the outer part In^ = J2ri>R IVnA'^ = 1 ^ will be 
due to the edge contribution. The radius R has been 
chosen to be about of one lattice constant smaller than 
the radius of the maximum circle that can be inscribed 
in a given GQD. Then, the state will be referred to as 
the edge state if In'' > In'' ■ Otherwise, we refer to it as 



to the bulk state. Wave functions of the typical edge and 
bulk states are illustrated in Fig. [1] 

The bulk and the edge states distinguished by the 
above criterium are also separated in energy. Namely, 
the edge states normally possess the energy \E\ < A, 
while the energy of bulk states \E\ > A. We will re- 
fer the energy interval of 2A around the Dirac point as 
pseudogap. It turns out that the value of pscudogap is 
approximately equal for all GQDs, characterized by the 
same inner radius R (see Fig. The pseudogap scales 
as A (X Jo/VN (X 7oa/i?. 

Two types of the edge states can be discerned: (i) the 
zero-energy states (ZES) that are degenerated and lo- 
cated exactly at E' = 0, i.e. in the middle of the pseudo- 
gap, and (ii) the dispersed edge states (DES) that have 
the non-zero energies, are symmetrically distributed with 
respect to E = 0, and fill the pseudogap. 

As it was shown by the graph-theory^ the total num- 
ber of ZES is related to the imbalance between the A and 
-B-typc atoms in the graphene flake: 

Vo>\Na-NbI (6) 

where the equality is taking place for the geometry of 
equilateral polygons. For hexagons 770 ~ 0, there are no 
ZES and all edge states are of the DES type. Contrary, 
for the equilateral triangles all the edge states are of ZES 
type and their degeneracy number is given by^ii 

= VnTs - 3. (7) 

Usually there are only few ZES for circular and randomly 
shaped GQDs. 

The number of ZES does not depend on magnetic held 
and therefore these levels do not contribute to the orbital 
part of susceptibility ([5]). In contrast, they are responsi- 
ble for the spin-provided super-paramagnetic response of 
ensemble of clusters in case of the half-fllled 7r-band when 
the Fermi energy is pinned at = 0. Indeed, according 
the Hund theorem, the number of single occupied states 
of degenerate level should be maximal, providing the to- 
tal uncompensated spin defined by the Lieb's ruleiiii^ 
S = ^ryo that brings the substantial contribution to the 
temperature-dependent spin-Curie paramagnetism. 

On the other hand the location of DES in hexago- 
nal, circular and random GQD depends on the applied 
field and therefore these levels are responsible for diamag- 
netism of graphene clusters, as it is calculated below. 

IV. HEXAGONS, CIRCLES AND RANDOM 
QUANTUM DOTS 

The susceptibility was calculated for zig-zag edge 
hexagons and circles of about of ten different sizes, hav- 
ing an inner radius in the range of 2 — 7 nm. For random 
quantum dots, the averaging has been performed over 
three different ensembles, characterized by mean inner 
radii Ri = 2.6 nm, R2 — 3.65 nm and R3 = 4.69 nm. The 
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FIG. 2: Electronic density of states and the diflerence of 
the edge and the bulk contributions intensities I^C^ — 7^''' (see 
text) as a function of energy for the hexagonal (a), circular 
(b), triangular (c) and random (d) GQDs. The shaded region 
indicates the pseudogap where the edge states are located. 
The levels have been artificially shown as Gaussian peaks with 
the dispersion 0.05eV. 



these geometries is vanishingly small that makes the dia- 
magnetic contribution dominant for all sizes. 



FIG. 1: Amplitude of the wave function for bulk (a-d) and 
edge (e-h) states in GQD of different shapes. We also show 
the maximal inscribed circles that have been used to discern 
between bulk and edge states. 



magnetic field varied between and 5 T, range in which 
the susceptibility remained approximately constant. All 
the plots are presented for i? = 5 T. Fig. [5] shows the 
density of states as a function of Fermi energy for hexag- 
onal and random-shape GQD. The shaded area indicates 
the region where the edge states are located. 

Fig. 13] shows the magnetic susceptibility per unit of 
area as function of Fermi energy for GQD of different 
sizes. As was qualitatively explained above, the diamag- 
netic peak of width 2A appears when the chemical po- 
tential crosses the pseudogap. This peak becomes wider 
with decreasing of GQD size. Beyond this zone, the or- 
bital susceptibility is a highly fluctuating function of the 
Fermi energy that oscillates between paramagnetic and 
diamagnetic sign. These oscillations have been recently 
interpreted for graphene ribbons, as a result of the sub- 
band structures^. 

As it was mentioned above, the number of ZES for 



V. TRIANGULAR QUANTUM DOTS 

All the edge states in triangular GQD with zig-zag 
edges are of ZES-type that do not change their energy 
as function of the field and do not contribute to the dia- 
magnetic susceptibility. The ZES for GQD of different 
sizes are located in the middle of the gap 2A as is shown 
in Fig. m We note that 2A indicates the real gap in the 
spectrum for the triangular GQD, and the value of the 
this gap is inversely proportional to the number of ZES 

y(FT3)-3' 

where the numerical constant is C ~ 5.56. 

The magnetic susceptibility of triangular GQD x^rb 
shown in Fig. [3] for nine inner sizes of 2.5 — 7nm. ft is 
provided by the out-of-gap delocalized electronic states 
and does not depend on ^ within the gap 2A because of 
absence of DES. These results match the analytical calcu- 
lations of Ref. [20||2^ for an infinite graphene sheet with 
the band gap 2 A, according to which the diamagnetic 
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FIG. 3: Susceptibility for triangular (a), random (b), circular 
(c) and hexagonal (d) GQDs of sizes 2 — 7nm. The shaded 
region indicates the (pseudo)gap. 



susceptibility per unit area is 
0(A-H) 



2A 



a 



? 9 2 



(9) 
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FIG. 4: The band gap for triangular GQD with different num- 
ber of atoms N . 



evaluated per unit area as 



Aspm 



(10) 



2 is the g-factor 



where /is is the Bohr magneton and g 
of the electrons with spin s = 1/2. 

In the opposite case of the strong Coulomb electron 
correlations, according to the Lieb theorem^, all 77^ ^ 
1 ZES electrons form the total spin of the cluster S = 
7yo/2. The super-paramagnetic susceptibility of ensemble 
of triangular GQD becomes even stronger, 



,^ _ 1 igfiBYSiS+l) 

Xspin _ 



(11) 



3 k^T ScrkBT 

The actual value of paramagnetic susceptibility should 
be somewhere in between of x^in ^^'^ Xs^in- 

Using Eqs. ([T5|) and (fTTj) we compare the spin- 
paramagnetic and orbital-diamagnetic contributions, 
presenting their ratio for the case of strongly correlated 
electrons as 



V-spm 



Xorb 



T.leV 



NksT 



(12) 



It follows that varying the size and temperature of the tri- 
angular quantum dots, we can expect the paramagnetic- 
diamagnetic crossover. In particular, for a temperature 
of T = 77 K, we predict that triangular quantum dots will 
be paramagnetic for the inner radius below i? ~ 97 nm. 



where 9(x) is the step function. 

Although ZES give no contribution to the orbital sus- 
ceptibility, they can be responsible for the huge paramag- 
netism provided by rjQ uncompensated electron spins lo- 
cated on the degenerate ZES levels. This happens in the 
case of small positive chemical potentials if the energy 
of electron-electron repulsion in each zero-energy state 
Ue-e > SO that thcsc levels remain half-filled. The 
corresponding Curie-type temperature-dependent para- 
magnetic susceptibility for non-interacting electrons is 



VI. SIZE DEPENDENCE 

We calculated the gap-zone-integral x{^)d^ of the 
shown in Fig. [3] diamagnetic susceptibilities for GQD of 
different shape and found that it preserves its intensity 
while the width 2A vanishes as the size of GQD increases. 
In the limit of infinite cluster with A — > this gives the 
McClurc (5-pcak of the graphene orbital susccptibilityi^: 

X{^i) = -aSifi). (13) 
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FIG. 5: Pseudogap A dependence on the number of atoms 
in GQD of different geometry. The inset shows the energy 
dependence of susceptibility for a circular GQD with A'' — 709 

The shown in FiglS] dependence of A on the number of 
atoms, is very similar for all shapes and follow the Eq.® 

The size-dependence of the orbital susceptibility for 
hexagonal, circular and random GQDs at /i = is shown 
in Figini It satisfies an empirical relation 

Xorb = PN^ (14) 

with A ~ 0.4 and /3 — 0.31a/7o. Although the orbital 
diamagnetism is originated from the edge currents of the 
low-energy DES, their perimeter contribution ^ N^/"^ 
can be reduced by the armchair or/and zigzag type 
boundary irregularities as well as by the wave function 
vanishing at the corners of GQD that explains the reduc- 
tion of the exponent index A slightly below 1/2. 




1000 2000 3000 4000 5000 
Number of atoms 

FIG. 6: Dirac-point orbital susceptibility at = 0. The solid 
lines correspond to the fit according Eqs. (|15|) and (|14|l . 



So far, there are no experimental measurements of 
magnetism of GQD but the susceptibility we obtained is 
in the range of the measured values for carbon nanotubes 
and buckyballs^. For convenience we also presented the 
susceptibilities in Figs. [3] and [6] in the experimental units 
of emu/g. 



VII. CONCLUSIONS 

Magnetism of GQD is provided by the edge states 
whose energy is located within the finite-size quantization 
pseudogap. The structure of the edge state spectrum and 
magnetic response of GQDs being strongly dependent on 
the geometric shape of the cluster. 

For hexagonal, circular and random GQD the edge 
states are dispersed within the pseudogap. Their posi- 
tion depends on the applied field, providing the substan- 
tial diamagnetic response of GQD. The diamagnetic sus- 
ceptibility as function of the chemical potential presents 
a peak of constant intensity, centered around ^ = 0. The 
maximum of the peak increases with GQD size whereas 
its width decreases, approaching the (5-function of Mc- 
Clure for infinite sheet of graphene. 

For triangular GQD the edge states are located exactly 
at the middle of the gap with the high degeneracy factor 
rjQ given by Eq.® that increases with size of the cluster. 
The zero-energy position of these levels do not change 
with the field and the diamagnetic response of triangular 
GQD ^OTb expected to be small. In a contrast, the 
uncompensated spins of electrons localized at ZES can 
provide the huge paramagnetic temperature-dependent 
contribution x^i„ of the Curie type. By comparison 
of susceptibilities Xmh ^'^d x^im Eq.(fT2|) we expect to 
have the crossover from paramagnetic to diamagnetic re- 
sponse in ensemble of triangular clusters with increasing 
the temperature and/or the GQD size. 

The strong dependence of magnetic properties of 
GQDs on their geometry, size and temperature provides 
the natural way to separate the graphene clusters accord- 
ing to their shape and size by application of the appro- 
priately designed non-uniform magnetic field and tem- 
perature cycle that can trap the different GQD in differ- 
ent points of space. It would be interesting also to study 
the specially cut nano-clusters of highly ordered pyrolytic 
graphite that can contain the separate graphene sheets 
with Dirac-like spectrum^ and therefore can have the 
similar magnetic behavior. 



For triangular GQD the size-dependence of the orbital 
susceptibility can be obtained by substituting the value 
of A from Eq.® into Eq.® at /i = 0: 



Xorb 



70 



2C 



(15) 



This estimation agrees with our numerical calculations 
as shown in FiglBl 
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